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ABSTRACT 


D diffuser throat diameter 


An interactive in viscid core flow-boundary layer 
method is presented for the calculation of turbomachine 
channel flows. For this method, a one-dimensional 
inviscid core flow is assumed. The end-wall and blade 
surface boundary layers are calculated using an inte- 
gral entrainment method. The boundary layers are 
assumed to be collateral and thus are two-dimensional. 
The boundary layer eguations are written in a stream- 
line coordinate system. The streamwise velocity pro- 
files are approximated by power law profiles. 
Compressibility is accounted for in the streamwise 
direction but not in the normal direction. 

Eguations are derived for the special cases of 
conical and two-dimensional rectangular diffusers. 

For these cases, the assumptions of a one-dimensional 
core flow and collateral boundary layers are valid. 
Results using the method are compared with experiment 
and good quantitative agreement is obtained. 

NOMENCLATURE 

Ah unblocked flow area = (1 - B) A q 

Ag geometric area 

A 1 -A 3 terms in velocity gradient equation 

a speed of sound 

B blockage 

b passage height 

Cf skin friction coefficient 


D 1 -D 2 terms in velocity gradient equation 
F entrainment function 

6 - 6 , 

H. shape factor, — 

1 Op 

<S ^ 

H 9 shape factor, — 

C op 

K 1 -K 3 terms in area equation 
K t , Ky curvature terms 
L diffuser length 

M relative Mach number, W/a 

M u wheel speed Mach number, U/a 

m exponent in velocity profile 

m mass flux 

P pressure 

PSTD standard pressure, 101 325 N/m 
Re Reynold's number 

Rg gas constant 

r radius 


specific heat at constant pressure 


"pr 


P - Pi 


pressure recovery coefficient, p , p - 

r le ~ 1 


s entropy 

T temperature 

T$td standard temperature, 288.2 K 
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t tangential coordinate 


INTRODUCTION 


U wheel speed 

W relative velocity 

x streamline coordinate 

y normal coordinate 

z blade number 

a angle between meridional streamline and 

axis of rotation 

$ angle between relative velocity vector and 

meridional plane 

y specific heat ratio 

5 boundary layer thickness 


6 1 

5 2 

e 

x 

M 

p 

T 

4> 


displacement thickness, / (1 


rt 


momentum thickness, 


f *L 

/ w 

1 e 


channel divergence angle 
prewhirl 

coefficient of viscosity 

density 

shear stress 



entropy function, 



dT 

T~ 


gj rotational speed 

Subscripts: 

A spanwise direction 

e edge of boundary layer 

r radial direction 

t,T tangential direction 

w wall 

x streamline direction 

y normal direction 

1 channel inlet 

2 channel exit 
Superscripts 

' absolute total conditions 

11 relative total conditions 


In the preliminary design of turbomachinery it 
is necessary to estimate the losses in the individual 
components. For rotors, stators and impellers, this 
is usually done by using empirical correlations for 
channel losses (e.g., skin friction, diffusion, etc. 
(!-£.))• For diffusers (conical, rectangular) this is 
done by empirically relating the diffuser pressure 
recovery coefficient to the diffuser geometry and 
inlet flow conditions (5,6>) . A more rigorous way of 
estimating these losses would be to calculate the 
surface boundary layers and from these obtain a mass- 
averaged estimate of channel losses. However, many of 
the present methods for calculating boundary layers 
(e.g., the code, BLAYER, ref. 7) are not ideally 
suited for calculating the boundary layers in turbo- 
machine flow channels. 

First of all, most boundary layer codes are 
written for external flows, thus, they require a 
specification of the free-stream velocity distribution. 
For internal flows, the velocity distribution is not 
known apriori and thus must be determined by an iter- 
ative procedure along with the continuity equation for 
the assumed inviscid core flow. This can be a very 
time consuming procedure especially for preliminary 
design calculations where several levels of iteration 
must be performed. 

Also, most boundary layer codes are written for 
two-dimensional flows along flat surfaces and thus 
ignore any coordinate curvature terms (i.e., metrics). 
This can lead to large errors for some internal chan- 
nel flows in which the channel end walls are converg- 
ing or diverging; e.g., nozzles and diffusers. 

There have been several methods proposed to cal- 
culate the performance of straight-walled, rectangular 
type diffusers (8-11) ♦ These methods are all two- 
dimensional in nature (i.e., they ignore end-wall 
effects) and are limited to the computation of incom- 
pressible, straight center-line diffuser flows. These 
methods are highly developed and predict the perform- 
ance of this class of diffusers very well even into 
the stalled regions of the diffuser maps. However, 
since all these methods are iterative, viscous- 
inviscid interaction methods, their running times can 
be much longer than that of the method proposed here. 
This could make them unsuitable for use in compressor 
preliminary design codes. 

For these reasons, a noninterati ve integral- 
entrainment method for calculating the flow fields in 
turbomachines has been developed. The method simul- 
taneously calculates the free-stream velocity distri- 
bution from an assumed one-dimensional core flow 
continuity equation together with the equations 
describing the surface boundary layers. Thus, only 
one pass through the flow channel is required thereby 
making the method very fast. The primary assumptions 
made for the boundary layer calculations are that the 
flow is basically two-dimensional and that the bound- 
ary layer profile can be described by a power law. 

The exponent in the power law is allowed to vary as 
prescribed by the equations cf motion. Compressibil- 
ity is accounted for in the streamline direction but 
not in the direction normal to the flow surfaces. 
Turbulent flow is assumed everywhere using the semi- 
empirical Ludwieg-Ti 1 Iman relation for skin friction. 
Approximate methods for calculating through regions of 
fully-developed or separated flows are also presented. 
Comparisons are made with experimental data for the 
cases of conical and rectangular diffusers and good 
agreement is demonstrated for cases where no separa- 
tion is indicated. Fair agreement is demonstrated in 
separated regions of the flow. 
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METHOD OF ANALYSIS 

The method of analysis used to calculate the flow 
surface boundary layers in the channel is an inteqral- 
entrainment method. The method is similar to that 
used for the calculation of end-wall boundary layers 
in Ref. 12. For this analysis, the core flow is 
assumed to be one-dimensional and is described by the 
one-dimensional continuity and momentum equations. 

Two-Dimensional Boundary Layer Method 

The following assumptions are made for this 
analysis : 

(1) The flow is steady and turbulent everywhere 

(2) The fluid is a perfect gas 

(3) The flow in the boundary layers is two- 
dimensional, i.e., there is no crossflow component of 
velocity. Therefore, the boundary layers are also 
collateral, i.e., the flow angle, 8, is constant 
through the boundary layer. 

(4) The channel flow surfaces are adiabatic 

(5) There is no density variation in the direc- 
tion normal to the channel boundaries although the 
density is allowed to vary in the streamline direc- 
tion, i.e., p = P (x). This approximation seems to 
work well for Mach numbers up to 1.0 and slightly 
beyond (12). For high supersonic flows, some other 
formulation would have to be used. 

(6) The boundary layer velocity profile in the 
streamline direction may be described by the power-law 
equation 



where W e , 5, and m are functions of distance along 
a streamline, x, only. The power law approximation 
adequately predicts pressure recovery up to separation, 
although the actual profile deviates from the power 
law approximation more significantly as separation is 
approached. 

To completely solve the problem, the boundary 
layer thickness, 6, and the boundary layer exponent, 
m, for each surface boundary layer calculated; the 
free-stream velocity, W e ; and the free-stream pres- 
sure, density and temperature must be determined. 

Thus, a total of (4 + 2n) equations must be solved 
where n is the number of surface boundary layers. 

If a streamline coordinate system is chosen (with 
x in the streamline direction, y in the normal to 
the wall direction, and t in the crossflow direction 
completing the orthogonal triad), the boundary layer 
equations for motion can be written as follows: 

Continuity 


9(pW ) a(pW ) 


ax 


3y x ' y 


pW x (K + K t ) = 0 


Momentum 


3 W 3 W 

PW X + pW y -jJT ~ P r “ 2 Sln 0 C0S 6 


= _ IP + 

3 X 


( 2 ) 


i 

r 3y 
(3) 



( 4 ) 


In Eq. (2), the curvature term, Ky, is defined as the 
curvature of the normal coordinate, y, in the 
t = constant plane and the curvature term, K^, is 
defined as the curvature of the tangential coordinate, 
t, in the y = constant plane. The normal coordinate 
is assumed to be everywhere normal to the flow. 

If Eqs. (2) and (3) are integrated through the 
boundary layer with the additional assumption that K t , 
Ky, a, 8 and r are independent of the normal 
coordinate, y; the following equations result: 


f s 3(pw x ) r s 

i ^ dr 'l 


3(pWJ 

— r — — dy + (K + 
3.y y 


V I •"« 
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3( P w; 


3 X 


— dy ♦ P W H + p(K + K.) I W‘ dy 


e ye y t 


f 


dy = 0 
(5) 


- prsu 2 sin a cos s = - s 4^- - t 

3X Xu) 


( 6 ) 


Equation (3) can also be evaluated for the isentropic 
core flow which results in the following equation 


d£ „ 

dx 


dW e 2 

eW e "3x + pra 


sin a cos 8 


(7) 


One additional unknown was added in the integration of 
Eq. (3), i.e., the normal velocity at the edqe of the 
boundary layer, W ye , as shown in Eq. (6). This can 
be modelled using Head’s entrainment relation as was 
demonstrated in Ref. 12. The resulting relation is 
gi ven by 

F 5 ( w e - W ye) = °-° 306 (H 1 " 3r °' 653 (8) 


Using the velocity profile from Eq. (1) and performing 
the indicated differentiations and integrations, 

Eqs. (5) to (8) can be combined to yield the following 
equations 


37=- 2m < 1 + 


m)(l + 2m) 


1 _^e . F 1 + 2m r 
W e dx Te 4ma L f 


(9) 


jds 

dx 


:i - 


M e’ 


dW 


dx 


1 

nn^ 


dm 


d* 




sin a cos 8 „ 2 , w 
r M u - (K y 




( 10 ) 


Equations (9) and (10) can be solved if the variation 
in free stream velocity, W e , is known. A detailed 
derivation of equations similar to Eqs. (9) and (10) 
is shown in the Appendix to Ref. 12. 

In order to solve these equations, the skin fric- 
tion coefficient, Cf, and the entrainment function, 

F, must be evaluated. The skin friction coefficient 
is obtained from the Ludwieg-Ti liman relation which is 
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-0.268 


C, - 0.246 Re. 
f 5 2 


-1.561 H 2 


(ID 


where H 2 = 1 + 2 m. 

The entrainment function is a function only of a 
boundary layer shape factor and was shown in Ref. 12 
to be gi ven by 


/ \ 0.653 

F = 0.0306 l-pr-i) (12) 

The variation in the free stream velocity, W^, 
can be determined if a one-dimensional core flow is 
assumed. In this case, the velocity can be determined 
from the continuity equation for the core flow which 
can be written 


m = pW e Ap (13) 

where A 5 represents the unblocked channel area 
normal to the flow. 

Using the results from Eq. (7) and differentiat- 
ing Eq. (13) results in the equation 


1 dW e 1 ( 

W e (1 - Mg) \ 


sin a cos B 


1 dA b) 


(14) 


The static temperature of the core flow can be 
determined from the energy equation for the core flow 
which is 


the definition in Eq. ( 8 )), and the boundary layer 
thickness equals half the channel height, Eqs. (9) 
and (14) can be solved as before. 

It would also be useful to be able to continue 
the calculation beyond an indicated point of boundary 
layer separation. This point is usually determined bv 
a limiting value of the shape factor, H? = 1 + 2 m. 
Once again, the boundary layer Eqs. ( 2 ) to (4) become 
invalid when the boundary layer separates. Also, the 
assumed power law profiles cannot accurately describe 
a separated flow. However, calculations show that the 
pressure recovery and losses can be reasonably pre- 
dicted using some simplifying assumptions. These are 
as follows: 

(1) The boundary layer Eqs. (2) to (4) and Dower 
law profile remain valid for separated flow. 

(2) The skin friction coefficient qoes to zero. 

(3) The entrainment function reaches a limiting 
value chosen to be the value at the separation point. 

(4) The shape factor is allowed to increase 
beyond the separation point. 

With these assumptions, the calculations can 
proceed beyond points of separation or fully developed 
flow. 

Loss Calculations 

Since the flow is assumed to be two-dimensional, 
the losses can be determined by a simple mass averag- 
ing procedure. Since the channel exit static pressure 
is determined in the solution of the boundary layer 
flow, all that is needed is a static enthalpy corre- 
sponding to a mass-averaged flow in the channel in 
order to calculate the mass-averaged entropy rise from 
the equation 


T = T" 



where 


T" 


2 

x , 2u\ - (u>r ) " 

'l 2 Ep 


(IE) 


(16) 


The static pressure and density can then be determined 
from either Eq. (7) and the equation of state for the 
core flow, or from the isentropic relations for the 
core flow. These equations together with Eq. (14), 
and Eqs. (9) and (10) written for each boundary layer, 
then provide the (4 + 2n) equations needed to solve 
the channel flow field. 


Fully Developed and Separated Flow 

When the boundary layers from opposite walls of 
a channel merge, an isentropic core flow no longer 
exists and the flow is said to be fully developed. In 
this case, the thin shear layer assumption used in the 
derivation of Eqs. (2) to (4) is no longer valid. 
However, it would be useful to carry out the calcula- 
tions after the flow becomes fully developed. In 
order to do this, some simplifying assumptions must 
be made. These are as follows: 

(1) The boundary layer Eqs. (2) to (4) describe 
approximately the fully developed flow field. 

(2) The centerline of the flow remains isentropic. 
Although this is consistent with PrandtVs mixing 
length theory of turbulence; in reality, it is not 
valid due to a nonzero shear stress qradient at the 
center 1 ine. 

With these assumptions and the fact that the 
entrainment function equals zero (as can be seen from 
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APPLICATIONS OF THE METHOD 

The method has been applied to the cases of 
conical and rectangular diffusers. In these cases, 
the assumptions of a one-dimensional core flow and 
collateral boundary layers are valid at least for 
nonseparated flows. 

Conical Diffuser 

The basic geometry of a conical diffuser is shown 
in Fig. 1. In order to solve Eqs. (9), (10) and (14), 
the cross-sectional area, Ap, and the curvatures, 

K y and K t , must be determined. The flow area is 
assumed to lie on a spherical segment as shown which 
is approximately everywhere normal to the flow. The 
unblocked flow area is then given by 



and the area term in Eq. (14) becomes 
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( 20 ) 


K 



1 “ (b - v« 1A ) 


The curvature terms for this case are qiven by 


Then for this case, since M u = 0, and 


K - 1 

2 (ex - a 1T ) 

k - £ 

j - (ex' - « 1T ) 


d{ l _ m ds + 6 dm 

“ 1 + "> (i + m) 2 ' 3 * 


( 21 ) 


The curvature terms for the tangential boundary layers 
are given by 


Eq. (14) becomes 



K t = 


(25) 


Id dx ” 0, 
e 1 


(22) 


The curvature terms for the spanwise boundary layers 
are given by 


(2 + 3 K,6, ) ? o c f 

where = - 2K^m F + K^(l + 2m) y- 

D 1 58 ^ - M e^ X + K lV + 2 K i 6 l ( 1 + "X 1 + 

Equations (9) and (10) complete the system of 3 
equations in 3 unknowns (W e ,m,6). The equations are 
solved using a fourth-order Runge-Kutta technique (13) . 

Results using this method have been compared with 
experimental data from a wind-tunnel diffuser section 
(14) and with data for conical diffusers (M. 

Figures 2 to 4 show comparisons of velocity profiles, 
static pressure coefficients, and displacement thick- 
nesses at the exit of the cross-leq diffuser of 
Ref. 14 for three different inlet Reynold's numbers 
and Mach numbers. In all cases the agreement is qood. 
For these cases, the boundary layers remained attached 
according to the analysis. Figures 5(a) to (d) shows 
comparisons of pressure recovery coefficient for sev- 
eral arbitrary cases of inlet Mach number, blockage, 
and cone angle for the conical diffuser data of 
Ref. 5. Almost all of the data lie in a region where 
the flow would be expected to be separated according 
to the analysis. However, the agreement between 
experiment and analysis is still fair. 


Rectangular Diffuser 

The geometry for a rectangular diffuser is shown 
in Fig. 6. For this case, the boundary layers in the 
tangential (blade to blade) and spanwise (hub to 
shroud) directions cannot be assumed to be equal due 
to the fact that the initial boundary layer thick- 
nesses may be unequal and the flow may become fully 
developed in one direction only. Thus, the tangential 
and spanwise boundary layers must be calculated sepa- 
rately. However, the core flow is still assumed to 
be one-dimensional . 

The unblocked flow area for this case is qiven by 


A b = (b - 26 1A )(2ex - 2« 1T ) 
and the term in Eq. (12) becomes 


i_^b 
A b dx 


d5 lT 

K 3 ” K 2 “dx 


d 6. 


I A 


dx 


(23) 


(24) 


where 


K y = 0; K t=7 

Equation (14) can then be expressed 

1 dW e \ + A 2 + A 3 

W e + °2 

where Aj = -2K ? m^F T - 2K 1 m^F A 


(26) 


(27) 


A 2=4(l + 2m T) 2 c fT + 2 1 U-2m A 2 )c fA 


A - K K * 1T K 5lA 
A 3 - " K 3 " K 2 ~ ~ K 1 T~ 


D 1 = (1 " M e><l + VlT - K 1#1a ) 

°2 = 2K 2 S 1T (1 + "V)(l + 2"V) + 2K 1 « 1a (1 + m A )(l + 2m A ) 


Equations (9) and (10) written for both the tangential 
and spanwise boundary layers then complete the system 
of 5 equations in 5 unknowns (W e ,m/\, s^my , . 

Results using this method have been compared 
with experimental data for rectangular diffusers (_6) . 

A comparison of calculated and experimental pressuFe 
coefficients for several arbitrary cases of inlet Mach 
number, blockage, and divergence angle, e, is shown 
in Fig. 7(a) to (d). The cases shown correspond to an 
inlet aspect ratio of one; however, the method is 
applicable for any arbitrary inlet aspect ratio. The 
cases were calculated assuminq equal boundary layer 
thicknesses and exponents on the inlet vane surfaces 
and end walls. The results show fair agreement 
between calculation and experiment even though the 
calculations indicate that almost all of the data 
shown lie in a region where the flow would be 
separated. 

CONCLUSIONS 


An analytical method has been developed to cal- 
culate boundary layers and losses in turbomachine 
flow channels. Results using this method have been 
compared with experimental data for the special cases 
of conical and rectangular diffusers. Good anreement 
is demonstrated for the flow regimes where the anal- 
ysis indicates no separation. Fair agreement is 
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demonstrated for the pressure coefficients calculated 
in the separated regions of the flow. 

Since the method is a noninterati ve integral 
method, it is very fast compared to other methods of 
analysis. Thus, it can be easily adapted for use in 
preliminary design or off-design codes for turbo- 
machinery components. 
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FIGURE 1.- CONICAL DIFFUSER GEOMETRY. 


W/We 





FIGURE 2.- BOUNDARY LAYER VELOCITY PROFILES AT EXIT OF 
.1 SCALE CROSSLEG DIFFUSER. 
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FIGURE 6.- RECTANGULAR DIFFUSER GEOMtTRY. 
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(B) M 1 = .6, B 1 = .12. 

FIGURE 7.- PRESSURE RECOVERY COEFFICIENTS FOR VARIOUS 
RECTANGULAR DIFFUSERS. 
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(C) M 1 = 1.0, B 1 = .06. 



FIGURE 7.- CONCLUDED. 
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